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ABSTRACT 

In the two decades since the first extra-solar planet was discovered, the detection and characteri- 
zation of extra- solar planets has become one of the key endeavors in all of modern science. Recently 
direct detection techniques such as interferometry or coronography have received growing attention 
because they reveal the population of exoplanets inaccessible to Doppler or transit techniques, and 
moreover they allow the faint signal from the planet itself to be investigated. Next-generation stellar 
interferometers are increasingly incorporating photonic technologies due to the increase in fidelity of 
the data generated. Here, we report the design, construction and commissioning of a new high contrast 
imager; the integrated pupil-remapping interferometer; an instrument we expect will find application 
in the detection of young faint companions in the nearest star- forming regions. The laboratory char- 
acterisation of the instrument demonstrated high visibility fringes on all interferometer baselines in 
addition to stable closure phase signals. We also report the first successful on- sky experiments with 
the prototype instrument at the 3.9-m Anglo- Australian Telescope. Performance metrics recovered 
were consistent with ideal device behaviour after accounting for expected levels of decoherence and 
signal loss from the uncompensated seeing. The prospect of complete Fourier-coverage coupled with 
the current performance metrics means that this photonically-enhanced instrument is well positioned 
to contribute to the science of high contrast companions. 

Subject headings: Instrumentation: interferometers, high angular resolution, techniques: high angular 
resolution 



1. INTRODUCTION 

The overwhelming majority of the 800 or so extra- 
solar planets confirmed to date*" have been detected us- 
ing the well-established transit ( [Burocki et a l.| |2010|) 
or Doppler (Tinney et al. 2002]) techniques. Although 
these indirect detection techniques are prolific in detec- 
tion yield, both are strongly biased to detecting large 
planets in close proximity to their parent star. Such 
selection biases have delivered a highly skewed picture 
of the mass/size/orbit distribution of extra-solar plan- 
ets which limits our understanding of planetar y evolu- 



tion and lea ves formation models unconstrained (Marois 
et al. 2010"). Unlike most indirect methods, high con- 



trast imaging techniques encompassing both coronogra- 
phy fGuyon 2003^ and interferometry ( [Angel fc Woolf] 
19971 ,Tuthin et al]|2000| make direct use of the light 
from the companion itself for scientific investigation. As 
a result exoplanetary companions can be studied in any 
orbit/period (Serabyn, Mawet & Burruss 2010) down 



Electronic address: jQvanovic.nem@gmail.com| 

^ MQ Photonics Research Centre, Dept. of Physics and Astron- 
omy, Macquarie University, NSW 2109, Australia 

Macquarie University Research Centre in Astronomy, Astro- 
physics and Astrophotonics, Dept. of Physics and Astronomy, 
Macquarie University, NSW 2109, Australia 

Australian Astronomical Observatory (AAO), PO Box 296, Ep- 
ping NSW 1710, Austraha 

^ Sydney Institute for Astronomy (SIfA), Institute for Photonics 
and Optical Science (IPOS), School of Physics, 
University of Sydney, NSW 2006, Australia 

^ Centre for Ultrahigh Bandwidth Devices for Optical Systems 
(CUDOS), Sydney, NSW, Austraha 

^ Observatoire de Paris, 5 place Jules Janssen, Meudon, France 

^ http://exoplanet.eu 



to and even below an angular separation corresponding 
to the diffraction limit (~ \/D where A is the wave- 
lengt h of light and D the diameter of the telescope mir- 
ror) ( Huelamo et al.|[2QTT ). In addition, it is possible to 
study planetary systems with arbitrary orbital inclina- 
tion to the observer, giving access to a larger population 
of potential targets. However, it should be noted that 
these techniques are currently limited to observations of 
young massive and/or actively accreting objects (warm 
and bright) while indirect methods can be used to detect 
older planets (cooler and hence fainter). 

One p owerful interferome tric technique is aperture 
masking ( Tuthill et al.|[2QQ0 ) which offers precise cal- 
ibration of wavefront structure for angular separations 
from the target star of around \/D. This is achieved 
by recording Fizeau interferograms generated by a non- 
redundant sparse-aperture mask placed in a re-imaged 
telescope pupil-plane. Such a scheme has been shown 
to be highly robust against the degrading effects of 
atmospheric phase-noise, particularly when the self- 
calibrating closure phase o b servable ([Ba ldwin 1986) is 
utilised (Lloyd et al. 20061 [ Lacour et~a lj |2011D . This 



technique has been responsible for a se r ies of recent high 
profile discoveries ( jKraus et al. 2008 2011[ ), including 
the first two newborn sub stellar companions captured in 
the pro cess of formation ( Huelamo et al.| [2QTT |Kraus fc" 
"Ireland 2012). Such young systems are valuable as they 
offer insights into planetary formation and allow us to 
constrain planetary evolution models. Despite these suc- 
cesses, aperture masking has several limitations includ- 
ing low throughputs (typically ~ 5-10% of the stellar 
photons are transmitted by a mask), incomplete Fourier 
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coverage, leakage of residual atmospheric noise due to 
phase corrugations across each sub-aperture and the two- 
dimensional nature of the interferograms makes poly- 
chromatic setups involving cross-dispersion difficult to 
realise optically. 

2. HARNESSING PHOTONICS 

Photonic technologies promise dramatic gains in the 
precision with which an interferometer can be fabricated, 
the fidelity of the data recovered, as well as in the ro- 
bust, stable performance over a wide range of challeng- 
ing environments. A recent photonic-reformulation of 
aperture masking replaces the mask with an optical sys- 
tem which feeds a number of single-mode optical fibres 
that sample the pupil-plane (Chang et al. 1998; Per- 
rin et al.|[2QQ6l|Lacour, Thiebaut fc Perrin 2007; TuthiH 



ing. A promising technology for achieving the desired 
3-dimensional wavegui de structures required is the laser 



et al. 2010); a scheme that addresses many of the cur- 
rent shortcomings. Firstly, because the light within a 
single-mode guide propagates with planar wavefronts, an 
interference pattern formed between the beams of two 
such guides will have a greater signal-to-noise ratio of 
fringe visibility (contrast) than if the atmospherically 
corrupted wavefronts were combined 100 x increase 



was demonstrated for the FLUOR instrument (Coude 



pu Foresto 1994)). Secondly, fibres can be routed from a 
2-dimensional pupil-plane to any configuration required 
at the output; an operation referred to as pupil remap- 
ping. This allows for the fibres to be reformatted for 
non-redundant beam recombination at the output while 
sampling the entire (highly redundant) telescope pupil 
at the input, enabling all available starlight to be used 
(up to 20 X increase in throughput compared to aper- 
ture masking if light is injected with a high-Strehl ratio 
adaptive optics (AO) system). A filled entrance pupil 
also yields complete Fourier sampling; a dramatic gain 
over a sparsely sampled mask. By reformatting waveg- 
uides into a linear output array, starlight ca n be injected 
into an integrated-photonic beam-combiner ()Berger et al 



2001 

cross-dispersed 



Laurent et "aL||2002[ [Benisty et 'al] |2009p , and/or 
lispersed for tne extraction of simultaneous spec- 

A key advantage of con- 



tral and spatial information 
ducting beam combination o n-chip is the possibility for 
nulling (Labadie et al. 2007), whereby the signal from 
the bright stellar source is cancelled out in one output 
port so that the faint signal from a companion can be 
detected (Bracewell 1978). 

The challenge with any optical stellar interferometer, 
including a pupil-remapper, is that the paths through 
the instrument must be matched to within the coher- 
ence length of the light {Lcoh. ^ 50 /im is typical for 
the bandwidths employed in the near-infrared). This is 
difficult to achieve with optical fibres because environ- 
mental infiuences such as temperature and strain must 
be carefully controlled over lengths up to several meters. 
Indeed, an optical fibre based pupil remapper known 
as the FIRST instrument has been demonstrated, but 
maintaining stable performance has p roved challenging 
dHuby et alT]|2012[ [Kotani et ar]|201QD . An ideal alter- 
native IS to use mtegrated-photomc components where 
all the waveguides can be embedded in a single minia- 
turized monolithic chip. In this way all the paths expe- 
rience a near-identical environment, with precision con- 
trol at fabrication allowing for tailoring of the trajectory 
of each guide in order to achieve path-length match- 



direct- write technique ( Nolte et ar][2003 Thomson et al 
[20TTI). In essence this involves sculpting tracks of modi- 
fied material, which act as waveguides, within the bulk of 
a dielectric (e.g. fused silica) wit h the use of a tightly fo - 
cused femtosecond-pulsed laser ( jGattass fc Mazur|2008 ). 

Here we present an integrated pupil-remapping inter- 
ferometer (IPRI) known as Dragonfiy, for high contrast 
imaging. This new instrument concept is ultimately in- 
tended to target ambitious scientific goals in exoplane- 
tary science. The instrument relies on the unique syn- 
ergy between several micro-scale optical technologies in- 
cluding the first photonic path-length-matched pupil- 
remapper, a steerable segmented-mirror array and micro- 
lens arrays. This paper presents the conceptual and prac- 
tical design features of the instrument, and summarises 
the results from the laboratory tests and the successful 
commissioning at the Anglo- Australian Telescope in sec- 
tion [3j From these results we extrapolate a predicted on- 
sky performance for the instrument when used in tandem 
with an AO system in section [4] In section [5j we round- 
out this body of work by outlining the narrow but crit- 
ical parameter space at small angular separations from 
the parent star, where the instrument will have few com- 
petitors and be able to make a key contribution to the 
field of exoplanetary science. 

3. EXPERIMENTS 

An artist's impression of the IPRI is shown in Figure [l] 
Starlight from a celestial target is collected by the tele- 
scope and routed to the photonic pupil-remapping chip 
into which it is injected by a micro-lens array (hexagonal 
lattice, 30 jam pitch). Although there were no specific 
limits to the spectral content of the signal that was cou- 
pled into the remapper, the IPRI was optimised for op- 
eration across the H-band. Each micro-lens of the array 
had a one-to-one correspondence with each of the 37 el- 
ements of a segmented mirror upstream and waveguides 
of the chip downstream. Since the segmented mirror and 
the micro- lens array were positioned in the conjugate 
pupil planes of the telescope respectively, the tip/tilt 
functionality of each mirror segment was used to steer 
the sub-pupil beams carefully in order to fine-tune the 
coupling into each waveguide. The mirror could also be 
used to compensate for path-length mismatches via its 
piston functionality. Although this was not utilised in 
the on-sky observations, it proved a very valuable fea- 
ture for our laboratory testing campaign. To minimise 
the amount of stray light injected into the pupil remap- 
per a mask was placed in close proximity to the front of 
the segmented mirror (i.e. close to the pupil plane, not 
shown in Fig. [T]). The mask was tilted at an angle to 
reject light from the segments that were not used and 
had circular sub-apertures (90% of the size of a mirror 
segments) to transmit the light for the segments that 
were used. After passing through the waveguides, a sec- 
ond micro- lens array (hexagonal lattice, 250 fim pitch) 
was used to re-collimate the emergent light before it was 
cross-dispersed by means of a 60° Flint glass prism and 
focused with anamorphic optics onto a research grade 
InGaAs camera (20 jam pixels). The linear dispersion 
was 29 nm/pixel (at 1.55 /im) giving a resolving power 
{R = A/ A A where A A is the bandwidth of the light), per 
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Fig. 1. — Schematic diagram of the IPRI. Light from the celestial target (depicted as a star with an orbiting planet in this figure) is 
collected by the telescope before propagating through the instrument and finally being collected at the detector. Orange circles indicate the 
correspondence between the segments used on each optical component. The diagram is figurative only: the aspect ratio and scale are not 
preserved. To follow the beam path downstream we have: CT - celestial target, TEL - telescope, COL - collimating lens, RAM - right angled 
mirror, SM - segmented mirror (IrisAO - PTTlll), RLl - relay lens 1 (200 mm focal length achromat), RL2 - relay lens 2 (10 mm focal 
length aspheric), MLAl - micro- lens array 1 (30 /im pitch), PRC - pupil remapping chip, MLA2 - micro- lens array 2 (250 /im pitch), CLl 
- cylindrical lens 1 (200 mm focal length), CL2 - cylindrical lens 2 (20 mm focal length), P - prism, CAM - camera (Xenics-Xeva- 1.7-640). 
Inset: Photograph of the IPRI instrument on a 900 x 600 mm breadboard. 



pixel of 53. 

3.1. The Pupil- Remapper Prototype 
3.1.1. Design and Fabrication 

A prototype pupil-remapping chip was designed and 
is depicted in Figure 2] It consisted of 8 waveguides 
which were constructed by interpolating with cubic spline 
functions between nodal points placed along each path. 
By moving the nodal points it was possible to attain 
physical path length matching between all waveguides 
to within 100 nm while maintaining a minimum sep- 
aration between guides of 30 /im in order to minimise 
cross-talk (Charles et alJ2012j). This design process only 
matched the path lengths physically because the optical 
contributions, such as the reduced effective index that a 
mode experiences as it propagates around a bend, are 
minor contributors that did not increase the path-length 
mismatch beyond the coherence length of the light. The 
design proved to be extremely challenging and the re- 
sulting device, is to the best of our knowledge the first 
3D, path- length matched chip consisting of unique routes 
(i.e. each waveguide is unique and not simply a mirror 
of any other). 

The device was inscribed into a boro-aluminosilicate 
glass substrate (Corning Eagle2000) with dimensions of 




Output 



Input 



Fig. 2. — A CAD rendering of the 3D paths taken by the 8- 
waveguide pupil remapper. Note: the aspect ratio for the drawing 
greatly inflates the transverse coordinates to exaggerate the bends 
and curves in the guides. Insets show images of two of the waveg- 
uides illuminated with red laser light. 

30 X 20 X 1.1 mm (L x W x i^) by the laser direct- 
write technique. An ultrafast titanium sapphire oscil- 
lator (Femtolasers GmbH, FEMTOSOURCE XL 500, 
800 nm centre wavelength, < 50 fs pulse duration) with 
a 5.1 MHz repetition rate, was used for inscription. The 
laser was focused into the sample using a 100 x oil im- 
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mersion objective lens (Zeiss N-Achroplan, NA= 1.25, 
working distance of 450 /im) . In order to prevent the im- 
mersion oil boiling during writing, the average power of 
the laser was controlled by using an external electro-optic 
pulse picker (BME KG) to reduce the repetition rate to 
1.28 MHz. Pulse energies of 160 nJ were used in conjunc- 
tion with a translation velocity of 8.3 mm/s in order to 
create waveguides that supported a single-guided mode 
at 1.55 jam. A set of Aerotech, air-bearing translation 
stages were used to smoothly translate the sample in 3 
dimensions. The entire device was written within half 
a minute. The device was then ground and polished to 
reveal the waveguide ends. As a result of the high rep- 
etition rates used all waveguides present ed here in were 
written in the cu mulative heating regime ( Eaton, Zhang] 
' fc Herman ||2QQ5l). 

The 8 single-mode waveguides sampled the re-imaged 
2D conjugate pupil plane of the telescope which had a 
diameter of 210 /im at the Coude focus (see Figure [3|. 
The starlight was reformatted into a ID equidistantly 
spaced (250 jam) output array that spanned a total of 
1.75 mm. In Figure |2] it can be seen that a lateral side 
step was used. This was done in order to prevent un- 
guided stray light reaching the output face where the 
waveguides termin ate which would ind uce phase errors in 
the measurement (iNorris et al. 2011). The distribution 
of the waveguides at the input oi the remapper (Figure [3| 
was chosen to conform to the hexagonal lattice of the el- 
ements of both the segmented mirror and the micro-lens 
array (overlaid in Figure Isl). A projection of the pupil 
of the primary mirror has also been overlaid in Figure [3j 
The inscribed circle of each hexagonal segmented mirror 
element and its corresponding micro-lens, subtended an 
area with 56 cm diameter on the primary mirror. For 
this first prototype a subset of 8 of the available lenslets 
was used in order to facilitate the characterization of the 
instrument, with the specific pattern chosen to provide 
adequate Fourier sampling. However, utilizing all avail- 
able lenslets is possible (with some work) and will be the 
focus of future prototypes and indeed the final instru- 
ment. 

3.1.2. Characterisation 

The prototype pupil-remapper was designed to oper- 
ate around 1.55 fim in the astronomical H-band. The 
slightly elliptical waveguides (13 x 14 ± 1 jam) are de- 
picted in Figure [3] and consisted of non-step index pro- 
files with peak index contrasts of > 6 x 10~^. Before 
the waveguides could be used for broadband interfero- 
metric applications, it was important to first determine 
the extent of the single-mode regime which is where spa- 
tial filtering could be exploited, as described above. The 
single-mode cutoff of the waveg uides was determi ned by 
the transmitted power method (Lang et al. 1994). This 
involves exciting the modes of the structure and inves- 
tigating the broadband transmission characteristics to 
identify the point in the spectrum where the losses drop 
abruptly (i.e. the transmission recovers). This point 
signifies that a second mode is now bound to the guide 
and can carry energy with relatively low loss which at 
longer wavelengths it was radiating away. The light from 
a broadband (super continuum) light source was injected 
and collected from the waveguides by butt coupling the 
device to single-mode optical fibres (SMF-28). It should 




Fig. 3. — Cross-sectional micrograph of the input facet of the 
pupil-remapping chip. Eight waveguides are clearly visible each 
of which consists of a complex index profile made up of a high 
index central region (bright central spot, positive index contrast > 
6 X 10~^) bordered by a low index region (dark spot, negative index 
contrast of < — 1 x 10""^) surrounded by a high index ring (weak 
positive index contrast < 3 x 10""^). The optical field at 1.55 /im 
is guided over the entire index profile. The image is overlaid with 
a hexagonal grid depicting the overlap with the elements of the 
segmented mirror and the injection micro-lens array. The dashed 
red overlay shows the projection of the3.9 m primary telescope 
mirror (the pupil). 

be noted that the fibres were not aligned for optimal 
transmission but rather they were offset laterally by 6 /im 
in order to excite the higher order mode of the guide. The 
transmitted spectrum was recorded by an optical spec- 
trum analyzer (OS A). The spectrum was normalised to 
that which was transmitted across a butt couple between 
the probe and collection fibres without the chip in be- 
tween. A normalised transmission spectrum for one of 
the waveguides from the remapper is shown in Figure [4] 
with the cutoff wavelength highlighted as the point where 
there was a dramatic recovery in the transmission. The 
mean single-mode cutoff wavelength for the 8 waveguides 
was 0.95 ± 0.03 jam which enabled the IPRI to operate 
in a highly broadband fashion across the entire Y, J and 
H bands (0.96-1.80 /im). It should be made clear that 
the high losses shown in Figure [4] are not indicative of 
the waveguides performance but are the result of the off- 
axis injection of the probe light in order to excite the 
higher order modes required for the measurement of the 
single-mode cutoff. 

The throughput of the pupil-remapper was charac- 
terised using the same fibre optic circuit used for the 
single- mode cutoff measurements, with a few slight mod- 
ifications; the light source was exchanged for a narrow 
linewidth laser diode operating at 1.55/im and a Ge- 
photodetector and power meter were used in place of 
the OSA. High precision, flexure translation stages were 
used to align the fibres with the waveguides, one at a 
time. The maximum transmitted power was recorded 
and normalised with respect to the power transmitted 
across a butt couple between the probe and collection 
fibres without the chip in between. Index matching im- 



Starlight Demonstration of Dragonfly 



5 




1000 1100 1200 
Wavelength (nm) 



1400 



Fig. 4. — Normalised broadband transmission spectrum for one 
of the waveguides in the pupil-remapping prototype. The single- 
mode cutoff wavelength is the point where the transmission recov- 
ers abruptly which has been highlighted in the Figure. The sharp 
features in the spectrum between 1.05 and 1.10 /im are due to the 
pump laser and the associated Raman lines of the super continuum 
source which were not entirely normalised. 

mersion oil was used at all interfaces to remove Fresnel 
reflections. The extent of cross-talk between the waveg- 
uides was examined by fixing the probe fibre for maxi- 
mum alignment into a particular guide while the collec- 
tion fibre was scanned between each of the waveguides 
at the output. The throughputs ranged from 5-50% for 
the 8 waveguides. The components of the loss are sum- 
marised in Table (T] It can be seen that ^ 20% is due 
to the non-negligible absorption by the Eagle2000 sub- 
strate ([Jovanovic et al. 2012 ) , 5% is due to a mode mis- 
match between the fibre and the waveguide modes (which 
was determined by calculating the overlap integral of the 
two fields), 30-93% is attributed to bend and transition 
losses, and < 0.1% is the result of cross-coupling. It is 
important to point out that the coupling loss reported 
was for the case when optical fibres were used to inject 
into the pupil-remapper, which is not the optical con- 
figuration for the IPRI. To account for this, the mean 
coupling efficiency was established with the micro-lens 
array to be 60 ± 2%. It should be noted that when 
the pupil-remapper was used in conjunction with micro- 
lenses, there was an additional Fresnel reflection loss at 
each facet of the pupil-remapping chip as ant ir effect ion 
(AR) coatings were not employed in this body of work. 
The higher than anticipated bend losses are the result of 
both the overly tight bend radii that were used in the de- 
sign of the prototype and the subsequent manufacturing 
errors that this caused when writing at high speeds. 

Although this level of performance was adequate for 
initial on-sky tests presented below, we believe we can 
minimise the losses of the pupil-remapper in future. This 
can be achieved by increasing the index contrast of the 
guides in order to optimise the coupling with the micro- 
lens array (80%), and minimise the cross-talk (< 0.05%). 
It should be made clear that the predicted maximum 
level of coupling of 80% is based on the limit when we 
consider the overlap between the Airy-shape d focal spot 
and the Gaus sian mode of the waveguide (Shaklan & 
Roddier 1988). In addition by utilising circular arcs in- 
stead of cubic spline functions and subsequently max- 
imising the bend radii, it should be possible to reduce 
the bend losses (< 15%) provided that waveguides are 
written at lower translation speeds to overcome manu- 



TABLE 1 

Breakdown of losses in the current pupil-remapping 
prototype and a prediction of the performance of an 

OPTIMISED DEVICE. ThE COUPLING LOSSES SHOWN ARE THE SAME 

FOR ALL waveguides/lenses, THE BEND/tRANSITION LOSSES 
REPRESENT THE RANGE OF THOSE FOR THE 8 WAVEGUIDES WHILE 
THE CROSS-COUPLING LOSS SHOWS THE MAXIMUM VALUE MEASURED 
AMONGST THE 8 GUIDES. ThE TOTAL THROUGHPUT OF THE 
CURRENT PUPIL-REMAPPER AS USED WITHIN THE IPRI IS 
SUMMARISED IN THE LAST LINE OF THE TABLE ALONG WITH A 
FORECAST FOR AN OPTIMISED SYSTEM. 



Component of loss 


Currently 


Predicted minimum 




measured 


achievable 




losses (%) 


losses (%) 


Absorption by substrate 


20 ± 1 


< 0.5 


Coupling losses 






Using optical fibres 


5±1 




Using micro-lens array 


40 ±2 


20 


Bend/transition losses 


30-93 


< 15 


Cross-coupling 


< 0.1 


< 0.05 


Total throughput of 






remapper within IPRI 


3-31 


> 68 



facturing difficulties. By using a highly transparent sub- 
strate like fused silica for example, the absorption losses 
in the substrate can also be reduced (< 0.5%) and finally 
by implementing AR coatings Fresnel reflections can be 
eliminated as well. With optimisation the throughputs 
of each sub-aperture of the rampper are projected to be 
> 68%. If we assume all other non- metallic optics within 
the IPRI are AR coated as well, we project a total sub- 
aperture throughput through the entire IPRI of > 50% 
should be achievable. 

3.2. The Integrated Pupil- Remapping Interferometer 

Prior to commissioning at the telescope, the pupil- 
remapper was aligned within the IPRI and tested in a 
laboratory environment with a 45 nm bandwidth light 
source. The laboratory tests were conducted without the 
cross-dispersing prism and anamorphic optics. Initially, 
the 7 possible interferometer baseline lengths were tested 
in turn, by misaligning the unwanted elements of the seg- 
mented mirror leaving only a pair of waveguides illumi- 
nated in each case. High visibility fringes were observed 
on all interferometer baselines, two of which are shown 
in the insets to Figure |5|a) and (b). The fringe pattern 
was windowed using a Hanning window and then inte- 
grated in the vertical direction in order to collapse each 
2-dimensional image into a 1-dimensional line profile. A 
Fourier transform was calculated for the line profile, the 
square of which yielded the power spectrum and squared 
visibilities (V^) of the pattern as shown in Figure [sj a) 
and (b). The average amplitude fiuctuation of the visi- 
bilities of the 7 interferometer baselines as a function of 
time was 2.2% RMS. 

Because the waveguides were equally spaced at the 
output of the pupil-remapper, a result of the fact that 
the remapper was designed to be used in tandem with 
a beam combining chip which takes this input format, 
a non-redundant subset was selected by the simple pro- 
cedure of misaligning unwanted channels with the seg- 
mented mirror array. An interferometer with 6 unique 
baselines was realised in this fashion using 4 waveguides, 
of which the associated fringe pattern, power spectrum 
and squared visibilities are shown in Figure [sjc). The 
theoretical maximum squared visibility is proportional 
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Fig. 5. — Power spectra displaying the detected Fourier compo- 
nents for the corresponding set of waveguides that were used (high- 
hghted by yellow dots while the unused waveguides are shown in 
blue). The insets show the recorded interference pattern in each 
case. The numbered parentheses indicate the baseline order with 
the shortest unit baseline corresponding to 250 /im (the separation 
between neighbouring waveguides). 

to where is the number of waveguides used 

(the basic math ematical f ounda tions of interferomet ry 

(20031)); 



Owing to its resilience against atmospheric and instru- 
mental phase noise, the closure phase has been the key 
observable for the detection of high c ontrast companion s 
with conventional aperture masking ([ Lacour et al.|20ir ). 
The closure phase is defined as the sum of the phase mea- 
sured around a closing triangle of baselines in the spatial 
frequency domain. A closing triangle comes about be- 
cause each aperture of the interferometer is responsible 
for 2 baselines and hence spatial frequency components 
and can be formed by using three sub-apertures, for ex- 
ample. It is insensitive to phase errors on individual 
channels, thus delivering an observable that is a function 
of the target object's intensity distribution only (Baldwin 
19861 . The stability of the closure phase as a function 



of time is of great importance as it is one of the param- 
eters (in addition to the Fourier coverage and others) 
that de fine the maximum c ontrast of a detectable com- 
panion ( [Lacour et al.| 2011). Stabilities at or below the 
1° level will see the IPRi reach scientific utility, while 
closure phases < 0.5° place the instrument on par with 
the present state-of-the-art achieved with aperture mask- 
ing in conjunctio n with AO systems ( Lacour et al.||2011[ 
Kraus & Irelandj|2012). The IPRI' s closure phase stabil- 



Lawson 



(2000) and Monnier 



ity was tested with artificially induced phase offsets (via 
pistoning of the segmented mirrors elements). It should 
be made clear that the pistoning was carried out with 
a single unidirectional sweep of the segmented mirror's 
elements, across their entire range, over a timescale of 
half-a-minute. This did not simulate the bidirectional, 
fast time scale phase fiuctuations induced by the atmo- 
sphere. A single frame was collected for each position of 
the mirrors' elements during the sweep. The visibilities 
were calculated as outlined previously. The closure phase 
was calculated for each frame by multiplying the complex 
visibilities of 2 of the Fourier components by the com- 
plex conjugate of the 3^^ that formed the closing triangle 
and then taking the argument of this value. The RMS 
fiuctuation in the closure phase was 0.4° over the pis- 
ton scan. As the segments were pistoned one-at-a-time 
over a range which corresponded to several multiples of 
the wavelength of the light, then the 0.4° RMS closure 
phase stability was measured with effective wavefront er- 
rors > 27r. However, current AO systems can deliver 
residual wavefront errors of ~ 250 n m while extreme AO 
systems aim for ^ 80 nm (SPHERE (Beuzit et aLjpOS ) 
performance fo r median seeing priv ately communicated 
by Boccaletti) ( Bouchez et al.|20Q9| ). Approximating the 
scaling of phase error s to be lin ear in the case of well- 



corrected wavefronts (Martinache 



are described in 

here this is V'^ = 6.063 (assuming even illuminationj and closure phase stabilities of 0.18° RmS and 0.058° RMS 



2010) yields expected 



is marked by a dashed line in Figure [sjc). It can be seen 
that the instrument generated high visibility fringes on 
all 6 baselines simultaneously which approach this value. 
The departure from the theoretical value for longer base- 
lines can be attributed to aberrations in the ad-hoc opti- 
cal setup downstream of the pupil remapper used to per- 
form this particular test which caused the higher fringe 
frequencies to wash out, and was not a problem for our 
on-sky testing which used the far superior optical setup 
shown in Fig. [l] Nonetheless the high visibilities ob- 
served confirm that the arms of the instrument were 
matched to significantly better than the coherence length 
of the light (A^/AA ~ 54 fim for the 45 nm bandwidth 
light source used). 



when used in tandem with an AO and an extreme AO 
system respectively (where we have also accounted for 
the effect of a random error on every sub-aperture si- 
multaneously). The closure phase stability projected for 
use with an AO system is a factor of 3 better than the 
0.5° threshold at which the IPRI becomes scientifically 
competitive while for extreme AO systems it is almost an 
order of magnitude below the state-of-the-art achieved in 
masking. We do not claim that this simple analysis can 
directly translate to the predicted on-sky performance 
where several sources of additional systematic error may 
arise. However this promising level of stability suggests 
that the instrument has the potential to compete with 
aperture masking as will be discussed in greater detail in 
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the later sections. 

3.3. On-sky observations 

The IPRI was tested on-sky on the nights of the 20^^- 
21^* of May, 2011 on the 3.9-m Anglo-Austrahan Tele- 
scope (AAT) at Siding Spring Observatory. As this 
telescope was not equipped with AO, the injection into 
the waveguides was severely degraded and interference 
fringes were in constant motion due to the turbulent 
phase perturbations imposed by the atmosphere. The 
median seeing for Siding Springs in May was 1.8 arcsec at 
0.55 jiia (from aao.gov.au) which resulted in a Fried pa- 
rameter, ro = 6.3 cm and more relevantly tq = 22 cm at 
1.55 /im.This value of ro was 2.5 x smaller than the pro- 
jection of each micro-lens on the primary mirror (56 cm) 
and hence the atmospherically induced wavefront errors 
were > tt radians across each micro- lens. This resulted 
in poor coupling of the starlight into the waveguides. 

In order to minimise the loss of wavefront coherence, 
short integration times were necessary but were limited 
by the readout noise of the InGaAs camera used (200 
electrons while a scientifically competitive HAWAII-2RG 
is < 18 electrons). As a compromise our observations 
were taken with 200 ms exposures, which were many 
times longer than the typical atmospheric coherence time 
(ro ~ 30 ms at 1.55 /im for Siding Springs). The target 
object was the red supergiant Ant ares, one of the bright- 
est objects in the near-infrared sky which was also at low 
air mass. A data set of 1000 consecutive frames was col- 
lected. An example single exposure (after cleaning and 
fiat-fielding) extracted from the data series is depicted 
in Figure |6] Clear fringes corresponding to the strongest 
baseline can be seen, with the slight tilt from vertical 
exhibited by the fringes with wavelength implying that 
the path-lengths through the interferometer generating 
that particular fringe were temporarily unbalanced (due 
to the varying atmospheric column depths above each 
sub-pupil). 

The cumulative power spectrum shown in Figure [7| 
was obtained by summing the squared one-dimensional 
Fourier transform over the 1000 frames. Two wavelength 
bands are clearly visible around 1.55 jam (H-band) and 
1.3 /im (J-band). Six peaks in the power spectrum can 
be seen in the H-band (marked by white circles), which 
show that fringes were recorded on all 6 baselines of the 
instrument. The system visibilities for these channels 
were calibrated by removing the effect of unequal beam 
fiuxes (which resulted from the unequal throughputs of 
the waveguides) in post processing, yielding a mean vis- 
ibility of 0.31. As no photometric signals were used to 
monitor the coupling into each waveguide, the rescaling 
was done based on the laboratory measurements of the 
throughputs, which was a reasonable first-order assump- 
tion. It should be made clear that the on-sky visibilities 
have been calibrated so that the maximum value each 
baseline can take is 1, which is in contrast to the labora- 
tory visibilities which were not calibrated. Up to 15% of 
the reduction in visibility can be attributed to the fact 
that Antares was partially resolved by the AAT (41.3 
mas angular size ( [Richichi fc LisilClQQOl )) leading to a re- 
calibrated mean system visibility of 0.36 for the longest 
baseline. 

A full numerical model of the telescope-atmosphere 
system was constructed to interpret the commissioning 
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Fig. 6. — A single image of a spectrally dispersed fringe pattern 
taken on-sky at the AAT while observing Antares. 
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Fig. 7. — Power spectrum of the interference fringes recorded by 
the instrument. The horizontal axis gives the spatial frequency, in 
inverse pixels, for the Fourier components. All fringes were well 
sampled. The white lines depict the theoretically expected spa- 
tial frequencies of fringes for our given remapping chip, reimaging 
optics and computed dispersion of the prism. The red and white 
circles highlight the detected spatial frequency components in J 
and H bands respectively. 

data. The system simulator included a single-layer phase 
screen with random fluctuations statistically matching a 
Kolmogorov turbulent spectrum. Adopting median Sid- 
ing Springs values for seeing and windspeed (7 m/sec 
aao.gov.au) in May, the code produced cross-dispersed 
interferograms. Light losses, beam asymmetries, pho- 
ton/readout noise, and the relatively long integration 
time of 200 ms were all incorporated in the model. Re- 
sults from the reduction of the simulated data gave me- 
dian system visibilities of 0.39 in the H-band, which were 
consistent with the on-sky value of 0.36 for the longest 
baseline. The visibilities had a signal-to- noise of about 
one from a single frame of data (both in simulations 
and as observed in the on-sky data), which is the the- 
oretically expect ed va lue when imaging in the speckle 
regime (Roddier) (19881) as was the case for our experi- 
ment due to the relatively long integration times. The 
simulator revealed that the majority of the degradation 
of the visibilities could also be attributed to temporal 
effects, where fringes were blurred due to the exposure 
time being significantly longer than the coherence time. 

It can be seen that several spatial frequency peaks were 
also recovered in the J-band (marked by red circles). Al- 
though our instrument was not specifically designed for 
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operation here, the short single- mode cutoff of the waveg- 
uides and the recovery of these spatial frequencies sug- 
gests a straightforward extension to the J-band in the 
future. 

Using the same code it was also possible to compare the 
measured throughput with expectations, with the find- 
ing that detected counts were about a factor of 2-3 lower 
than predicted. Such losses could easily arise if the see- 
ing was a little worse than simulated, and/or from known 
issues with the microlens alignment with respect to the 
chip. Recovering expected throughputs is a key com- 
ponent in pushing the science towards a fainter class of 
scientifically compelling targets. Although these results 
indicate there is still some work to do, they are encour- 
aging enough to project that high efficiency operation 
should be possible with a more mature optical setup. 

Closure phases were extracted corresponding to all of 
the 4 possible closing-triangles and we discuss the stabil- 
ity of the one with the highest signal-to-noise of these. 
Each of the 5 spectral channels spanning the H-band 
yielded closure phases of ± 5° standard error (zero is 
the expected closure phase for a nearly unresolved cen- 
trosymmetric object such as Antares), which were con- 
sistent with numerical simulations (0 ± 4° standard er- 
ror) that included seeing. We therefore conclude that 
the on-sky closure phase performance can be entirely 
attributed to uncorrected seeing (temporal fringe blur- 
ring) and hence there were no unexpected or unexplained 
drops in performance at the telescope as compared to the 
laboratory. 

4. PREDICTED ON-SKY PERFORMANCE WITH 
ADAPTIVE OPTICS 

With the closure phase stability for the on-sky tests 
explained by the presence of seeing, it is reasonable to 
use the closure phase stability of the laboratory measure- 
ments as a proxy for predicting the on-sky performance of 
the instrument if used in conjunction with AO systems. 
It is important to understand that the following extrap- 
olation of on-sky performance when used in conjunction 
with AO does not take into account all the systematic er- 
rors that are present when operating on a telescope and 
is meant only as a guide to the ball park performance 
that may be achieved in an ideal case scenario where the 
systematics can all be calibrated out. 

We can arrive at a predicted maximum contrast of de- 
tecting a companion as follows. As outlined previously, 
the laboratory measurement of the closure phase stabil- 
ity can be rescaled to account for wavefront errors across 
all sub-apertures simultaneously around a closing trian- 
gle with the magnitude equal to that expected of the 
residuals from an AO or an ex treme AO system (^ 250 
and ~ 80 nm ([Bouchez et al. |2009, ) respectively). This 
yields a 0.18° RMS and 0.058" RMS closure phase ex- 
cursion, respectively. This level of closure phase stabil- 
ity can be converted into a limit for which a companion 
could be detected at an angular separation oflX/D from 
the stellar host (pa rameterized for a 7- aperture system 
Lacour et al. ( 2011[ )) to give a con- 
.a) for an AO system and ~ 7000 : 1 



this means that a companion which is 750 times fainter 
than the parent star could be detected with a signal-to- 
noise of 3.3 when used in conjunction with an AO sys- 
tem. These values compare favourably with the current 
state-of-the-art contrasts for aperture masking that have 



using equation 2 in 
trast of ^ 2200 : 1 ( 

(Icr) for an extreme AO system. When converted to con- 
trasts for typical confidence levels (3.3cr, 99.9%) used in 
aperture masking observations we get ~ 750 : 1 for AO 
and ~ 2100 : 1 for extreme AO systems. For clarity 



been achieved; 100-400 : 1 (Huelamo et al. 2011 



Kraus 



Ireland|2012 j), but cannot be directly conipared to the 



aperture masking results as they were captured in longer 
wavelength bands. The predicted levels of contrast would 
currently limit the H-band version of the instrument to 
exploring the Brown dwarf regime. However, by moving 
further into the mid-infrared (L'-band) where the stars 
are fainter and the young planets are therefore relatively 
brighter, Jupiter-like planets may potentially be within 
reach. 

In practice, extrapolations of performance from the 
laboratory to the telescope are fraught with difficulty. 
For a start our above prediction might be regarded as 
pessimistic, for it strictly only applies to a small num- 
ber of detected frames and random errors in the closure 
phase noise could be reduced by taking a longer data 
run (there is a limit to this as well). However, an im- 
portant source of systematic error not considered thus 
far is that encountered when slewing between a science 
target and a point-spread-function reference star, and 
there may be persistent-speckle or non-common-path er- 
rors that come into play which degrade the maximum 
contrast achievable (although it is worth noting that ex- 
treme AO systems are specifically engineered to avoid 
such issues). It is also important to bear in mind that the 
system commissioned represents a very simple testbed; a 
final scientific instrument would benefit from two major 
scalings in data quantity (and quality). For simplicity 
of analysis, the closure phase statistics presented so far 
strictly apply to a single wavelength channel. Critical 
analysis of the ensemble of wavelength channels across 
the H-band not only boosts the SNR with the square 
root of sample size, but also delivers differential phase 
data. Zimmerman et al. have shown a factor of 3-4 
improvement in SNR by exploiting such spectro-spatial 
data (Zimmermann et al. 2011). The second major scal- 
ing is associated with the expansion of the array with 
a next-generation remapper capable of utilizing a larger 
fraction of the full pupil. For example a scheme which 
recombined 4 sets of 7 sub-apertures in non-redundant 
configurations which would be relatively simple to imple- 
ment would yield 60 independent closure phase measure- 
ments; a factor of two gain in contrast over the numbers 
quoted above. As our ^ 2000 : 1 projection is already 
competitive with the state-of-the-art over this range of 
spatial scales, we will (conservatively) adopt this value in 
the discussion of scientific reach which follows, but note 
that there are substantial grounds for the belief that a 
next generation instrument could do significantly better. 

5. DEFINING THE SCIENCE DISCOVERY SPACE FOR AN 
INTEGRATED PUPIL REMAPPING INTERFEROMETER 

Three parameters are required to define the niche pa- 
rameter space for exoplanets with the IPRI: the angular 
separation (between the planet and the star), the bright- 
ness of the star and the contrast between the planet and 
the star. The real strength of the IPRI lies in being 
able to detect companions at spatial scales of half to a 
few X/D from the parent star with moderate contrast 
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2000 : 1 at 3.3cr as described above). In compari- 
son, the Gemini Planet Imager (GPI) instrument which 
makes use of an extreme AO system an d a Lyot corona- 
graph will target detections at > 3 A/I) ( [Macin tosh et al.| 
|2006p . There are several sophisticated coronagraph de- 
signs such as the phase-induced amplitude apodisation 
(PIAA) coronagraph being employed in the SCExAO in- 
strument which aim to reduce this inner blind s pot to 
^ l\/D ( |Lozi, Martinache fc and GuyonJ|2009[) . Al- 
though contrasts as high as 10^ : 1 have been demon- 
strated in a laboratory environment, wavefront calibra- 
tion at separations around 1\/D have proven extremely 
difficult thus far and contrasts have not exceeded 100 : 1 
in this region to date. This means that our IPRI of- 
fers highly competitive contrasts for planet detection in 
the 1\/D separation region for the immediate future 
against extreme AO/coronograph technologies. This is 
an extremely important region of the angular separation 
space when we consider the nearest star forming regions 
(Scorpius-Centaurus at 380 — 470 lyrs and Taurus- Ariga 
460 lyrs). The solar-like stars in these regions are be- 
lieved to harbour many young planets which are still 
warm and hence very luminous from the formation pro- 
cess making them easier to detect. An angular separa- 
tion of 1\/D for an 8-m telescope in the H-band cor- 
responds to 7 AU (the approximate distance of Jupiter 
from the sun) at the distance to these star forming re- 
gions. This represents a sought-after parameter space 
as it corresponds to solar system scales (^10 AU from 
the parent star) and hence this forms the basis for the 
demand for instruments addressing the regime of small 
angular separations from the star. 

At present, aperture masking would represent the chief 
competitor to the IPRI at small angular separations. Al- 
though the predicted contrast for the IPRI is superior to 
masking when used with AO, the current throughputs 
are lower on an aperture-by-aperture basis. This, as out- 
lined in the experimental section above, we believe we 
can address in future to realise a remapping interferom- 
eter with > 50% throughput across each sub-aperture. 
However, when using the IPRI in tandem with a next 
generation extreme AO system, it is not the throughput 
that will determine the limiting brightness of the host 
star (magnitude) that can be studied, it is the level of 
signal required for the wavefront sensor of the AO system 
to lock successfully. Once the AO system is locked, long 
integration times can be used to maximise the signal- 
to-noise ratio of the Fourier observables (only limited by 
readout noise) . In this way the IPRI will extend even the 
reach of aperture masking in the region of small angular 
separations from the host star. 

In future the PIAA coronagraph being developed at 
Subaru telescope will also become a strong competitor 
in this angular separation regime. Although it currently 
offers th roughputs > 50% for separations > 2. 2\/D from 
the star (| Lozi, Martinache fc and Guyon|2009l ) , with new 
PIAA lenses it is projected to push this inner working an- 
gle to 0.6A/D. As the only other next generation high 
contrast imaging instrument being currently developed 
to image as close in as the IPRI to the parent star (ignor- 
ing aperture masking for the moment), it is interesting to 
note that the two instruments share a similar level of cost 
and complexity; the IPRI has a custom made photonic 
chip, a segmented mirror and 2 micro- lens arrays while 



the PIAA has two sets of custom made lenses, a spider 
removal plate and a binary mask. This reinforces the 
fact that the IPRI is uniquely positioned in the field of 
high contrast imaging with the few competitors required 
to exploit similar levels of sophistication to compete. 

The pupil remapping technology pioneered here also 
represents a basic enabling platform for future photonic 
instrumentation. With it, light from a telescope pupil 
may be injected into a planar waveguide structure while 
preserving the inherent wavefront coherence. Harness- 
ing the myriad of devices now in use from the photon- 
ics and communications industry for the manipulation of 
light in integrated chips, elements such as splitters, cou- 
plers, gratings and delay lines may all be incorporated 
into future instruments. One example of the use of such 
elements might be in the construction of a nulling inter- 
ferometer, an instrument designed to null out the bright 
star with high suppression leaving behind the signal from 
the faint companion. This interferometric analogue of a 
coronagraph would yield still greater advances into high 
contrast imaging. 

6. CONCLUSION 

We present the first experimental results from an inte- 
grated pupil-remapping interferometer from both labo- 
ratory and on-sky testing. The powerful combination of 
micro-scale optical technologies, namely the segmented 
mirror, microlens arrays and the pupil-remapping chip, 
enables advanced photonic processing and control of 
starlight in a simple and elegant fashion. The instrument 
has demonstrated high visibility fringes and stable clo- 
sure phases (0.4° RMS for > 27r wavefront error) from the 
ffist accurately path-length-matched, 3D photonic device 
consisting of uniquely routed waveguides in a laboratory 
setting. The on-sky tests showed performance metrics 
(throughputs, visibility and closure phase recovery) con- 
sistent with those expected given wavefronts corrupted 
by uncorrected atmospheric turbulence. We believe that 
with the use of extreme AO on larger telescopes both 
high visibilities and high levels of closure phase stability 
will be achieved. For a modest 7-aperture system, the 
ball-park projection of the faint companion contrast at- 
tainable is ~ 2000 : 1 (3.3<j) for an angular separation 
down to 1\/ D from the parent star. A key feature of our 
architecture is the ease with which it may be dramatically 
scaled by populating the entire pupil plane with waveg- 
uides. An IPRI offering near-complete Fourier coverage 
and spectrally dispersed data, coupled with predicted 
performance levels of closure phase stability from an ex- 
treme AO equipped telescope would deliver unsurpassed 
performance over the most critical region of angular sep- 
aration: that corresponding to solar system scales in the 
nearest star- forming regions. The possibility of such ex- 
citing performance illustrates the immense promise for 
new-generation astrophotonic instrumentation to make 
unique contributions to observational astrophysics. 
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